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This  study  investigates  and  demonstrates  the  enhancement  in  thermoelectric  power  factor  for  n-type 
AIZnO/InAIZnO  multilayer  quantum  wells  as  compared  to  their  counterpart  bulk  films.  A  10-20% 
improvement  is  observed  for  operating  temperatures  <700  °C.  Fabricated  structures  are  composed  of 
50  periods,  with  targeted  individual  layer  thicknesses  of  10  nm.  The  best  performing  multilayer  shows 
an  electrical  resistivity  and  Seebeck  coefficient  of  1000  pV/K  and  at  700  °C.  In  addition,  a  theoretical 
relationship  is  derived  between  the  thermoelectric  performance  and  correlating  microstructure  that 
demonstrates  the  deterioration  of  electronic  transport  properties  at  increased  interface  roughness  levels. 
To  determine  the  microstructure  and  interface  roughness  of  the  films,  X-ray  and  spectroscopy  techniques 
were  used.  The  proposed  model  is  based  on  how  quantum  well-width  fluctuations  caused  by  increased 
interface  roughness  leads  to  the  localization  of  carriers  and  thus  a  decrease  in  electrical  conductivity. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Today,  most  of  the  energy  resources  used  in  industrial  processes 
are  discharged  as  waste  heat  into  the  environment.  Such  exhaust 
heat  accounts  for  approximately  60%  unused  energy.  Hence,  ther¬ 
moelectric  (TE)  energy  conversion,  which  converts  waste  heat  into 
electricity,  has  received  much  attention.  The  conversion  efficiency 
of  a  thermoelectric  material  is  determined  by  the  dimensionless 
figure  of  merit,  ZT  =  oe2T/pK,  where  T  is  the  absolute  temperature, 
at  is  the  Seebeck  coefficient,  and  p  and  k  are  the  electrical  resistivity 
and  thermal  conductivity,  respectively.  However,  the  three  physical 
parameters  comprising  Z  are  all  functions  of  the  carrier  concentra¬ 
tion,  n.  Taking  a  and  p  to  be  directly  related,  and  p  and  k  varying 
inversely  to  one  another  (the  Wiedemann-Franz  law),  the  improve¬ 
ment  in  ZT  is  not  an  easy  task.  Currently  available  thermoelectric 
materials  have  a  ZT  of  less  than  ~1  and  the  device  efficiency  of 
producing  electric  power  rarely  exceeds  5%.  This  performance  has 
limited  TE  generators  to  niche  applications  where  requirements 
for  remote  operation,  reliability,  no  moving  parts,  and  silent  opera¬ 
tion  have  outweighed  the  more  negative  aspects  of  high  cost  and 
low  conversion  efficiency.  It  is  estimated  that  with  a  ZT  >  2  and  a 
temperature  difference  of  200  K,  a  10%  heat-to-electric  conversion 
efficiency  can  be  obtained  which  would  make  this  technology  more 
feasible  for  mainstream  applications  [1,2]. 

Some  thermoelectric  materials  being  researched  today  are 
nanocomposite  alloys  of  conventionally  known  TE  materials  such 
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as  Bi2Te3,  PbTe  and  SiGe  doped  with  additional  impurities.  An 
example  is  La-doped  PbTe-Ag2Te  nanocomposite  which  uses  the 
Ag2Te  precipitates  to  enhance  phonon  scattering  while  La-doping 
allows  for  carrier  concentration  control  and  results  in  ZT  values 
of  >1.5  at  775  K  [3].  A  secondary  set  of  materials  include  skutteru- 
dites,  which  use  a  mechanism  of  impurity  void  filing,  which  have 
resulted  in  n-type  ZT  values  of  up  to  1.7  for  CoSb3  compounds 
[4].  However,  these  materials  are  commonly  toxic,  low  in  abun¬ 
dance  as  natural  resources,  and  thus  not  environmentally  benign. 
Based  on  these  disadvantages,  metal  oxides  have  emerged  as  an 
alternate  potential  candidate  for  high  temperature  thermoelectric 
applications.  In  addition  to  their  temperature  stability  and  resis¬ 
tance  to  oxidation,  they  can  also  be  deposited  via  a  variety  of 
feasible  techniques  that  are  easily  scalable.  However,  the  main 
concern  is  that  conventional  metal  oxide  films  have  exhibited  rel¬ 
atively  low  levels  of  ZT  (~0.34  at  1000  °C)  due  to  poor  carrier 
mobility  from  their  ionic  structures  [5]. 

The  well  known  and  highly  investigated  p-type  superconductor, 
NaCo204,  was  the  first  metal  oxide  observed  to  show  promising 
thermoelectric  properties.  Upon  further  research,  the  layered 
structure  of  Ca3Co409  has  emerged  as  the  leading  material  candi¬ 
date  with  a  ZT  of  1. 2-2.7  at  600  °C  [6,7].  This  is  attributed  to  the 
alternate  stacking  of  hexagonal  Ca2Co03  layers  and  rectangular 
rock  salt  layers,  which  results  in  increased  phonon  scattering  and 
significantly  reduced  k.  A  second  approach  is  to  reduce  the  ionicity 
of  the  metal-oxygen  bonds  in  order  to  increase  the  carrier  mobil¬ 
ity,  an  example  being  ZnO-based  thermoelectrics.  Currently,  the 
most  widely  studied  n-type  material  is  Al-doped  ZnO  (AZO),  which 
has  obtained  a  ZT  value  of  0.1  at  1200  °C  [8].  More  recently,  by 
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means  of  co-doping  with  additional  low  level  impurities  (Ga  or  In), 
an  enhancement  of  ZT  was  found.  This  is  either  due  to  an  amplifi¬ 
cation  in  the  phonon  scattering  at  the  increased  number  of  inter¬ 
faces,  or  the  enhancement  in  the  electrical  transport  properties 
from  raising  the  density  of  the  states  which  coincides  with  the 
addition  of  carrier  electrons  [9,10].  However,  ZT  values  for  n-type 
films  remain  inferior  for  practical  applications. 

The  use  of  lower  dimensionality  materials  has  been  theoretically 
proven  to  be  an  underlying  mechanism  in  improving  the  thermo¬ 
electric  properties  of  materials,  predicting  ZT  values  >3  [11].  Some 
low  dimensional  structures  include  quantum  wells  (QWs),  superlat¬ 
tices,  quantum  wires,  and  quantum  dots.  Within  these  low-dimen¬ 
sional  regimes  where  quantum  effects  are  dominant,  the  energy 
spectra  of  electrons  and  phonons  can  be  controlled  by  altering  the 
size  of  the  structures  which  allows  for  the  manipulation  of  the  mate¬ 
rial  properties.  In  theory  this  regime  gives  rise  to  new  materials  even 
though  the  material  is  made  up  of  the  same  atomic  structure  as  its 
parent  bulk  material.  Multilayer  QWs  also  known  as  superlattices, 
consist  of  more  than  50  periods  with  individual  layer  thicknesses 
of  <20  nm  [12].  The  first  material  acts  as  the  active  layer  while  the 
second  material  with  a  larger  band  gap  acts  as  a  barrier  for  the 
charge  carriers.  The  band-gap  offset  thus  creates  a  well  that  is 
periodically  repeated  throughout  the  multilayer  [13].  Fig.  1,  is  the 
structural  design  of  our  samples  which  utilize  the  in-plane  transport 
of  heat  and  current  flow  (parallel  to  the  substrate),  whereas  conven¬ 
tional  TE-based  superlattices  highlight  the  reduction  of  thermal  con¬ 
ductivity  via  cross-plane  transport  flow  [14].  From  the  in-plane 
transport  methodology  we  benefit  from  (a)  enhanced  density  of 
states  near  Ff,  which  results  in  an  enhanced  Seebeck  coefficient 
and  (b)  an  increase  in  the  energy-dependence  of  the  carrier  mobility 
dju(E)  =  d E  which  results  from  increasing  the  scattering  time’s 
energy  dependence,  by  preferentially  scattering  electrons  depend¬ 
ing  on  their  energy  (quantized  energy  band  levels).  In  addition,  by 
using  the  in-plane  transport  structural  geometry  we  can  obtain 
much  greater  distances  between  heat  source  and  sink  which  allows 
for  a  lower  heat  flux,  qx  =  kAT/x,  where  AT  is  the  temperature  differ¬ 
ential  and  x  is  the  distance. 

To  obtain  feasible  devices  it  would  be  essential  to  scale-up  this 
technology  to  large  area  substrates.  For  this  reason  an  appropriate 
deposition  technique  must  be  chosen  and  verified.  Molecular  beam 
epitaxy  and  atomic  layer  deposition  are  deposition  techniques 
which  provide  excellent  precision  for  ultra  thin  layers  thicknesses, 
but  have  low  throughput  due  to  low  growth  rates.  Out  of  all  the 
techniques  available,  RF  sputtering  provides  high  quality  films  at 
high  growth  rates,  and  is  relatively  easy  to  scale  up  to  provide 
feasible  large  area  depositions  as  needed  for  alternative  energy 
applications.  However,  one  well  known  issue  coinciding  with  the 
high  growth  rate  of  RF  sputtering  is  the  increased  levels  of  surface 
roughness.  For  our  device  design  optimization  it  is  essential  to 
optimize  the  TE  properties  of  the  lOnm  QWs  which  means  we 
must  understand  the  impact  of  well-width  fluctuations  caused  by 
interface  roughness  (IFR).  We  have  developed  a  model  which  cor¬ 
relates  the  dependence  of  electrical  resistivity  to  IFR.  For  the  mul¬ 
tilayer  stacks  under  investigation  the  achievable  roughness  level 
from  deposition  becomes  extra  critical  because  initial  layer  rough¬ 
ness  will  subsequently  impact  the  layer  roughness  deposited 
thereafter. 


.  __ -  „  Multilayer  Thin  Films 

Heat  and  Current  Flow  In-Plane  ^ 

""^hot  a  —  "^"cold 

Glass  Substrate/ 

Fig.  1.  A  representative  diagram  of  our  multilayer  QW  thermoelectric  design,  using 
in-plane  heat  and  current  transport  (Not  to  scale). 


In  this  work  we  deposit,  characterize,  and  measure  the  thermo¬ 
electric  performances  of  multilayer  QWs  consisting  of  50  alternat¬ 
ing  periods  of  AZO  and  (In,  Al)  co-doped  ZnO  with  2-10  at.% 
indium.  Within  this  paper  we  will  refer  to  these  samples  as  InAZ02, 
InAZ05,  InAZ08,  and  InAZOlO  which  correspond  to  2  at.%  In  doped 
Al0.o2ZnO,  5  at.%  In  doped  Al0.02ZnO,  8  at.%  In  doped  Al0.02ZnO,  and 
10  at.%  In  doped  Al0.o2ZnO  films,  respectively.  In  addition,  we  will 
compare  thermoelectric  performance  to  previously  investigated 
bulk  films  of  the  same  materials.  Furthermore,  we  will  discuss 
our  newly  developed  theoretical  model  relating  the  electrical  resis¬ 
tivity  dependence  on  interface  roughness  and  the  subsequent  fit  to 
our  experimental  data  obtained  via  X-ray  analysis. 

2.  Experimental  details 

Bulk  and  multilayer  films  have  been  deposited  using  a  stainless  steel  physical 
vapor  deposition  reactor  equipped  with  two  confocal  magnetron  sputtering  guns. 
The  guns  were  driven  independently  by  13.56  MHz  rf  power  supplies.  A  computer 
controlled  automated  system  was  used  to  operate  the  opening  and  closing  of  the 
target  shutters  in  order  to  control  the  co-doping  process.  The  targets  were  commer¬ 
cially  sintered  3"  disks  of  A1203  (2  wt%)-doped  ZnO  (99.99%  purity,  Plasmaterials 
Inc.)  and  indium  (99.95%  purity,  Plasmaterials,  Inc.).  Preceding  all  depositions,  the 
reactor  was  evacuated  using  a  combination  of  a  mechanical  pump  and  turbomolec- 
ular  pump  to  a  base  pressure  of  <5  x  10-7  Torr.  In  order  to  vary  the  Indium  in  film 
concentration  between  2  and  10  at.%,  the  In  target  RF  power  was  increased  from  10 
to  20  W  accordingly  [14].  Preliminary  baseline  experiments  were  performed  on 
AZO  to  optimize  the  deposition  in  terms  of  high  crystallinity,  low  resistivity,  and 
low  root  mean  squared  (RMS)  surface  roughness.  The  process  parameters  are  Ar/ 
02  gas  ratio  of  3,  2  mT  operating  pressure,  100  W  AZO  target  power,  and  a  150  °C 
substrate  temperature  which  is  located  40  cm  from  the  target.  For  comparison  pur¬ 
poses  all  films,  both  bulk  and  multilayer  films  have  a  total  targeted  thickness  of 
1  pm. 

The  50  periods  of  AZO/InAZO  multilayer  QW  films  are  deposited  with  a  targeted 
individual  layer  thickness  of  10  nm,  because  this  will  result  in  quantum  well  effects 
(<20nm).  Deposition  times  varied  due  to  the  increased  growth  rate  (1.8-2. 7  nm/ 
min)  for  higher  indium  RF  power  (10-20  W)  [14].  All  films  were  deposited  on 
100  mQ-cm  n-type  (1 00)  Si  wafers  as  blanket  films  for  characterization  purposes. 
These  devices  were  fabricated  by  forming  1.5  x  0.5  cm  plateaus  of  multilayers  via 
reactive  ion  etching.  Chemical  selectivity  to  the  Si  substrate  was  not  an  issue  and 
did  not  result  in  any  substrate  pitting.  Devices  were  then  inserted  into  a  custom  de¬ 
signed  clamp  which  allows  for  the  Cr/Au  electrical  contact  deposition  on  the  short 
0.5  cm  ends  (Fig.  1)  without  any  subsequent  deposition  on  the  top/bottom  of  the 
sample.  This  was  done  on  the  samples  with  multilayer  structures  as  well  as  just 
on  the  substrate,  in  order  to  allow  for  the  separation  of  Seebeck  coefficient  measure¬ 
ments  [15].  This  will  be  explained  in  further  detail  in  the  Section  3.2.  In  addition,  the 
substrates  were  cleaned  using  the  standard  RCA  method,  and  subsequently  dried 
with  N2. 

3.  Results 

3.1.  Structural  and  morphological  properties 

Samples  deposited  consisted  of  50  periods  of  10  nm  layers  with 
QWs  ranging  from  0  to  10  at.%  indium,  resulting  in  a  total  thickness 
of  1  pm.  Phase  identification  was  carried  out  on  a  Scintag  X-ray 
diffractometer  (XRD)  equipped  with  a  Cu  Kal  X-ray  source  and  a 
horizontal  wide-angle  four-axis  goniometer  with  stepping  motors 
which  allowed  for  independent  or  coupled  6/29  axis  motion.  The 
multilayer  films  were  initially  annealed  for  2  h  at  900  °C  in  atmo¬ 
sphere  to  allow  for  necessary  re-crystallization  and  then  scanned 
from  20  to  60°.  The  collected  XRD  patters  were  compared  to  refer¬ 
ence  patterns  from  the  standard  Joint  Committee  for  Power  Diffrac¬ 
tion  Standards  powder  diffraction  file.  For  comparison  purposes,  all 
sample  data  has  been  normalized  to  its  maximum  peak  intensity. 

After  normalization,  the  results  obtained  are  identical  to  the 
bulk  XRD  patterns  obtained  in  a  previous  investigation  of  (In,  Al) 
co-doped  ZnO  bulk  films  [7].  The  XRD  pattern  in  Fig.  2,  show  that 
0  and  2  at.%  In  films  have  a  sole  detectable  diffraction  peak,  which 
is  nearly  consistent  with  that  of  the  (002)  standard  ZnO  crystal 
plane.  Using  the  Debye-Scherrer  formula  the  grain  sizes  for  the  0 
and  2  at.%  indium  films  are  32  ±  2  and  38  ±  2  nm,  respectively.  This 
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Fig.  2.  Bragg  diffraction  pattern  for  Al.02ZnO/InxAl.o2ZnO  multilayer  quantum  wells,  with  0  ^  x  ^  0.1.  Note  the  shift  in  preferred  film  orientation  from  ZnO  (002)  to  ZnO  (1 00) 
for  AZO/InAZ08. 


is  attributed  to  the  larger  atomic  radius  of  indium  being  incorpo¬ 
rated  into  the  ZnO  lattice  structure.  At  5  at.%  indium,  we  observe 
a  (1 00)  ZnO  peak  that  with  even  further  doping  becomes  the  pre¬ 
ferred  film  growth  orientation.  Also  at  8  and  10  at.%  indium,  In- 
based  peaks  materialize  which  is  a  sign  of  elemental  agglomeration 
within  the  crystal  structure. 

In  thin  films  analysis’s  the  surface  morphology  majorly  impacts 
the  use  of  spectroscopy  techniques,  and  for  this  reason  similar  we 
have  verified  similarity  using  atomic  force  microscopy  (AFM).  A 
Veeco  Dimension  3100  Scanning  Probe  Microscope  at  CNSE  was 
used,  along  with  a  sharp  silicon  probe  that  is  ideal  for  tapping 
mode  at  42  N/m  and  320  kHz.  The  average  surface  roughness 
(Ra)  was  measured  from  five,  1  pm  x  1  pm  size  scans  for  each 
sample  and  the  results  were  1.1,  1.3,  1.8,  and  3.1  nm  for  InAZ02, 
InAZ05,  InAZ08,  and  InAZOlO,  respectively.  For  all  images  taken 
(Fig.  3),  the  surfaces  appear  to  be  contaminant  free  and  without 
any  significant  peaks  or  valleys.  We  therefore  safely  assume  that 
surface  variation  will  not  impact  our  X-ray  results.  Also  apparent 
is  the  fact  that  the  AZO/InAZ02  appears  to  have  the  smallest  grains 
on  the  surface  which  is  in  agreement  with  the  XRD  data  obtained 
from  the  bulk  films.  This  may  be  attributed  to  the  lowest  lattice 
mismatch  between  the  barrier  and  well  layers  for  all  the  samples. 

In  order  to  characterize  the  microstructure  of  the  QW  multilay¬ 
ers  in  terms  of  periodicity,  layer  density,  and  interface  roughness, 
X-ray  reflectivity  (XRR)  profiles  were  recorded  by  a  Rigaku  Ultima 
III  X-ray  diffractometer  equipped  with  parallel  beam  optics.  Paral¬ 
lel  beam  optics  are  capable  of  delivering  high-intensity  X-ray 
beams  with  low  divergence.  The  X-ray  mirror  provided  a  high- 
intensity  beam  with  0.05°  divergence;  the  divergence  slit  was  set 
at  0.05  mm  to  minimize  the  sample  size  illumination  effects.  The 
acquired  XRR  profiles  were  analyzed  by  the  Motofit  software  pack¬ 
age.  The  sample  was  modeled  as  a  constant-density  multilayer 
stack  of  AZO  and  In-doped  AZO.  The  layers’  thickness,  density, 
surface  and  substrate-film  roughness  were  used  as  fitting  parame¬ 
ters.  The  software  performed  the  least-squares  regression  to  ex¬ 
tract  the  roughness  and  density  information  from  the  XRR  traces. 

In  Fig.  4,  we  observe  that  QW  samples  with  2-8  at.%  indium 
doping  within  the  well  demonstrate  an  ordered  layer  structure 
with  excellent  periodicity.  This  is  inferred  from  the  well-defined 


first  and  second-order  reflections.  However,  the  sample  with 
10  at.%  indium  doping  within  the  well  shows  evidence  of  non-peri¬ 
odicity  and  non-superlattice  behavior.  This  can  be  attributed  to 
elemental  agglomeration  of  the  indium  resulting  in  a  low  layer 
thickness  to  interface  roughness  ratio. 

In  order  to  verify  the  XRR  analysis,  cross-sectional  transmission 
electron  microscopy  TEM  microimages  were  taken  with  a  JOEL 
201 0-F  FEG  TEM  operated  at  200  kV,  shows  periodic  stacking  of 
the  constituent  layers  of  the  multilayer  film  as  shown  in  Fig.  5a. 
A  more  high  resolution  image  is  seen  is  Fig.  5b,  where  we  can  iden¬ 
tify  the  indium  doped  wells,  which  appear  as  dark  layers  due  to 
their  increase  in  density.  These  images  are  of  the  AZO/InAZ08 
sample  and  confirm  uniform  10  nm  layer  thicknesses. 

3.2.  Thermoelectric  performance 

In  order  to  calculate  the  thermoelectric  power  factor  (oe2/p)  of  the 
samples  both  the  electrical  resistivity  and  Seebeck  coefficients  are 
needed.  The  substrate  contribution  to  the  electrical  resistivity  is 
considered  negligible  due  to  the  resistivity  ratio  between  substrate 
and  film  being  >10000:1  [16].  The  resistivity  was  measured  both 
using  a  4-point  probe  at  room  temperature  and  a  2-point  probe  tech¬ 
nique  from  300  to  975  K.  Obtained  values  were  within  2%  of  each 
other  at  RT.  By  obtaining  the  total  resistance  measurement  with 
the  2-probe  technique  and  knowing  the  sample  geometry  we  can 
back-calculate  the  QW  structure  electrical  resistivity.  For  our  de¬ 
vices  the  measured  total  resistance  includes  the  QW  structure-me¬ 
tal  contact  resistance  and  the  in-plane  QW  structure  resistance.  For 
this  reason,  structures  of  different  lengths  are  used  to  eliminate  the 
effect  of  the  contact  resistance  and  to  obtain  the  sample  resistivity. 
Obtained  results  indicate  that  the  contact  resistance  values  are 
<0.1  mO-cm,  which  accounts  for  <5%  for  the  lowest  resistivity  film. 
Bulk  films  measurements  provided  are  the  results  of  a  previous 
investigation  that  used  the  same  measurement  techniques  [7]. 

The  in-plane  Seebeck  coefficient  measurement  is  straightfor¬ 
ward:  building  up  a  temperature  difference  across  the  sample  and 
measuring  the  output  voltage,  however,  in  our  device  the  Seebeck 
coefficient  of  the  substrate  (asut>)  will  contribute  to  the  overall  mea¬ 
surement.  In  order  to  separate  out  the  sample,  a  measurement  is 
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Fig.  3.  AFM  scans  for  three  multilayer  stacks  with  varying  levels  of  indium  doped  wells.  Results  show  all  films  are  free  of  large  agglomerations  and  that  grain  size  increases 
with  increased  doping  levels. 


Fig.  4.  X-ray  reflectivity  profiles  of  50  periods  of  AZO/InAZO  with  various  indium 
doping  levels  from  2%  to  10%.  Plots  indicate  defined  reflections  for  the  2,  5,  and 
8  at.%  indium  films  which  are  indicative  of  good  periodicity  of  the  multilayer 
structures. 

performed  on  both  the  substrate  as  well  as  the  combination  of  sub¬ 
strate  and  film  stack  (acombo)-  We  then  apply  Eq.  (1 )  to  our  results  in 
order  to  determine  the  Seebeck  coefficients 

^comboAT  —  ^MlATml  +  $subATsub  (1) 

The  Seebeck  coefficient  setup  includes  two  cylindrical  heaters 
that  are  bored  into  copper  blocks,  each  with  an  independent  tem¬ 
perature  controller  (Omega  Inc.).  Then  by  using  thermocouples 


placed  at  the  end  of  the  samples  a  50  K  temperature  gradient 
was  maintained,  with  the  hot  side  temperature  ranging  from  325 
to  975  K.  A  voltmeter  was  soldered  onto  a  nearby  position  which 
measured  the  voltage  differential  [17]. 

Fig.  6a-c  exhibit  the  temperature  dependent  measurements  of 
the  in-plane  electrical  resistivity,  Seebeck  coefficient,  and  thermo¬ 
power,  respectively.  Fig.  6a  and  b  are  experimentally  measured, 
and  c  is  a  calculation  based  of  the  prior.  Fig.  6a  has  been  plotted 
on  a  log  scale  in  order  to  show  the  full  range  of  measured  values. 
We  observe  that  p  for  the  QW  multilayer  structures  has  been 
reduced  by  up  to  an  order  of  magnitude  at  specific  temperatures 
(500-750  K)  and  doping  levels  (1  and  2  at.%  indium).  The  amount 
of  improvement  over  the  parent  bulk  films  varies  for  each  of  the 
QW  structures.  For  example  at  575  K,  the  electrical  resistivity  of 
the  QW  structures  as  compared  to  the  In-doped  AZO  bulk  film  is 
improved  by  63%,  89%,  91%,  and  50%  for  the  2,  5,  8,  and  10  at.% 
indium  doped  films,  respectively. 

Seebeck  coefficients  were  also  measured  at  50  K  intervals  be¬ 
tween  325  and  975  K  and  are  shown  in  Fig.  6b.  As  expected  all  sam¬ 
ples  have  negative  Seebeck  coefficients  demonstrating  n-type 
behavior.  The  bulk  and  multilayer  samples  both  show  a  linearly 
decreasing  a,  with  increasing  operating  temperature.  Unlike  the  p 
measurements  we  do  not  see  an  improvement  in  a  that  surpasses 
the  values  of  the  bulk  components.  The  resultant  a  value  is  actually 
an  intermediate  between  that  of  the  two  parent  films  it  is  composed 
of.  In  addition,  for  all  films,  we  see  that  with  increasing  at.%  indium 
concentration  there  is  a  trend  of  decreasing  a  values.  For  instance, 
at  575  I<  the  Seebeck  coefficients  for  InAZ02,  InAZ08,  and  InAZOlO 
are  -143,  -202,  and  -210  pV/K,  respectively. 

Based  on  the  electrical  resistivity  and  Seebeck  coefficient  mea¬ 
surements,  the  thermoelectric  power  factor  was  then  calculated 
and  plotted  against  temperature  in  Fig.  6c.  We  observe  that  the 
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Fig.  5.  Cross-sectional  TEM  images  showing  (a)  full  sample  QW  stack  and  (b)  partial  stack  w/dark  layers  indicating  indium  doped  quantum  wells. 


multilayer  structures  with  2  and  5  at.%  In  doped  wells  show  the 
most  promising  power  factors.  Specifically  within  the  temperature 
range  of  550  to  750  K  we  see  an  increase  in  power  factor  by  more 
than  3x  as  compared  to  the  parent  bulk  materials.  The  AZ0/InAZ08 
multilayer  shows  an  intermediate  level  of  performance  as 
compared  to  its  bulk  counterparts,  while  AZO/InAZPIO  shows  no 
improvement. 


4.  Discussion 

Our  experimental  results  indicate  that  the  sample  structures 
have  significantly  reduced  electrical  resistivity  at  low  levels  of  In 
doping  within  the  QW  (2  and  5  at.%)  as  compared  to  both  the  par¬ 
ent  bulk  films,  which  is  attributed  to  two  main  factors.  First,  a  por¬ 
tion  of  the  improvement  is  material  dependent.  Thru  the  addition 
of  a  more  electrically  conductive  film  (InxAl0.o2ZnO)  we  gain  a 
moderate  decrease  in  the  total  multilayer  structure  resistivity.  This 
is  proven  by  the  AZO/InAZOlO  sample  which  shows  no  sign  of 
layer  periodicity  by  XRR,  yet  measures  a  resistivity  of  50  mQ-cm 
at  575  K,  which  is  an  intermediate  values  of  the  parent  films,  105 
and  1.7  mQ-cm  for  the  In0.iAZO  and  AZO,  respectively.  Secondly, 
when  the  resistivity  is  lower  than  that  of  either  of  the  bulk  films, 
this  is  attributed  to  enhanced  quantum  confinement  effects  within 
the  quantum  wells.  For  example,  AZ0/InAZ05  measures  a  resistiv¬ 
ity  0.93  mQ-cm  at  625  K,  which  corresponds  to  a  55%  and  90% 
improvement  over  the  bulk  AZO  and  InAZ05  films,  respectively. 
From  Fig.  6a  we  also  note  that  similarly  to  all  the  bulk  films,  the 
resistivity  of  the  multilayer  structures  also  increases  with  higher 
operating  temperature  due  to  increased  scattering.  However,  at 
operating  temperatures  >650  K  rate  of  rise  (Ap/AT)  is  much  more 
significant  for  the  QW  structures  as  compared  to  bulk  films.  This  in 
theory  can  be  the  result  of  an  increase  in  multilayer  hopping 
caused  by  the  charge  carriers  gaining  enough  thermal  energy  to 
cross  the  band-gap  offset  of  the  quantum  well  [12]. 

In  order  to  determine  the  optimum  doping  level  for  the  QW 
structure  it  is  important  to  calculate  the  thermoelectric  power 
factor  due  to  the  inter-dependence  between  electrical  resistivity 
and  Seebeck  coefficient.  The  relationship  between  %  doping  in 
the  QW  and  electrical  resistivity  in  not  a  simple  straight-forward 
relationship  as  we  have  just  discussed  and  Fig.  6a  exhibits.  How¬ 
ever,  the  Seebeck  coefficients,  from  Fig.  6b  show  an  almost  linearly 
trend  of  improving  a  for  increasing  at.%  indium.  This  is  attributed 
to  the  increase  in  polycrystallinity  of  the  ZnO  matrix,  which  results 
in  a  reduction  in  carrier  mobility  for  both  phonons  and  electrons. 
This  has  been  verified  thru  XRD  analysis  (Fig.  2).  As  a  result  by  tak¬ 
ing  into  account  an  intermediate  level  of  doping  and  obtaining  high 
quality  films,  the  TE  power  factors  achieved  signify  the  advantage 


of  utilizing  quantum  confinement  for  AZO  and  (In,  Al)  co-doped 
ZnO  multilayers  at  high  operating  temperatures. 


4.1.  Interface  roughness  effect 


In  order  to  provide  appropriate  targets  for  device  fabrication  we 
need  to  understand  the  impact  of  material  quality  and  understand 
the  need  for  repeatability  and  fabrication  control.  From  our  exper¬ 
imental  results  we  have  identified  a  significant  decline  in  perfor¬ 
mance  for  multilayer  stacks  with  high  levels  of  IFR.  This  agrees 
with  the  theory  that  within  a  QW  any  deviation  from  a  perfect 
interface  should  impact  the  electrical  conductivity,  at  that  at  low 
temperatures  IFR  is  the  limiting  mechanism  for  carrier  transport 
[18].  The  following  model  has  been  developed  based  on  the  fact  that 
roughness  at  the  interface  leads  to  fluctuations  in  the  QW  width, 
which  in  turn  modulates  the  confinement  energies  and  results  in 
a  fluctuating  potential  for  the  2D  motion  of  confined  carriers.  The 
fluctuation  in  band  energy  is  given  by  the  partial  derivative  of 
energy  E,  with  respect  to  well  width  L 


*<r>  -  W m 


(2) 


where,  E  is  the  confinement  energy  and  is  given  as  ft2 IT2 /2mzL2,  and 
mz  is  the  z-direction  hole  effective  mass,  h  is  the  reduced  Planck 
constant,  and  A(r)  is  the  interface  roughness  at  a  position  r=  (x,y). 
The  influence  of  IFR  on  the  mobility  of  2D  electrons  is  never  precise 
since  the  roughness  itself  is  not  straight  forward  to  model.  This 
roughness  is  characterized  by  a  z-height  A  and  a  correlation  length 
A,  and  assumes  a  Gaussian  like  form  because  of  the  always  chang¬ 
ing  interface  structure 


(A(r)A(r'))  =  A2  exp 


\r-rf\ 

S  ) 


(3) 


Taking  into  account  that  the  interaction  is  actually  the  average 
over  all  possible  configurations  of  the  IFR,  the  correlation  function 
is 


mm) 


(4) 


where,  the  variance  of  the  function,  with  the  energy  squared,  is 


(5) 
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Fig.  6.  Plot  of  electrical  resistivity  (a),  Seebeck  coefficient  (b),  and  thermopower  (c)  vs.  temperature  for  bulk  and  multilayer  films  with  varying  at.%  indium.  The  multilayer 
films  with  observed  quantum  wells  show  up  to  a  3x  improvement  in  thermopower. 
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Correlation  Length  (nm) 


z-height  Roughness  (nm) 


Fig.  7.  A  3D  plot  of  resistivity  in  relation  to  z-height  roughness  and  correlation  length  for  a  well-width  of  10  nm  at  room  temperature.  Resistivity  increases  exponentially  with 
z-height,  and  linearly  with  correlation  length. 


Next,  by  using  a  previous  derivation  we  can  calculate  the  elec¬ 
tron  mobility  [19], 


25/2  e  / mz\2Cn]J 2 


(6) 


and  include  it  in  our  material  resistivity  calculation,  where  n  is  the 
electron  concentration,  p  is  the  hole  concentration,  and  jin  and 
are  the  electron  and  hole  mobility,  respectively.  For  our  calculation 
we  consider  the  impact  of  holes  negligible,  due  to  high  n  doping. 


P  qGu„n  +  jUpP) 

To  model  our  samples  we  use  roughness  RMS  values  obtained 
via  XRR  analysis.  The  simulated  data  gives  us  both  well-width 
and  z-height  roughness  (the  peak  height  in  the  Fourier  transform). 
We  estimate  the  correlation  length  from  the  assumed  Gaussian 
form  of  the  IFR.  Previously,  it  has  been  reported  that  the  IFR  can 
be  neglected  in  quantum  wells  over  4.5  nm  in  width  [20],  but  this 
is  only  true  for  a  low  level  of  roughness  (<1  nm),  which  would 
require  a  deposition  technique  other  than  sputtering.  More  impor¬ 
tantly  we  believe  that  it  is  the  ratio  of  interface  roughness  to  quan¬ 
tum  well  width  that  is  most  important. 

From  the  3D  plot  in  Fig.  7,  as  well  as  in  Eq.  (6),  we  can  see  that  a 
change  in  z-height  roughness  has  an  exponential  affect  in  deter¬ 
mining  the  sample  resistivity,  whereas  the  correlation  length  im¬ 
pact  is  linear.  It  is  observed  that  with  A  values  below  2  nm  it  is 
possible  to  achieve  resistivity  on  the  order  of  mQ-cm  for  all  A. 
Above  2  nm,  localized  areas  with  varying  potentials  will  begin  to 
develop  due  to  well-width  fluctuations.  This  difference  in  potentials 
leads  to  carrier  localization,  and  thus  higher  resistivity. 

To  complement  our  theoretical  findings,  AZ0/InAZ05  multilayer 
samples  were  analyzed  with  XRR  to  obtain  a  set  that  are  similar  in 
every  aspect,  yet  with  varying  IFR.  The  differences  in  IFR  can  be 
attributed  to  the  location  as  respect  to  substrate  origin  during  depo¬ 
sition.  From  best  fit  XRR  simulations  we  chose  four  samples  with 
the  properties  shown  in  Table  1 .  In  order  to  estimate  the  correlation 
length,  we  use  a  A  measured  from  XRR,  and  a  calculated  A,  based 
on  a  previously  determined  Gaussian  relationship  [12].  Both  the 
measured  and  calculated  values  are  shown  here  in  Table  2.  The  per¬ 


cent  error  between  the  measured  and  calculated  values  is  within 
10%,  except  for  sample  1019_c  which  has  a  substantially  higher 
level  of  roughness.  From  this  we  can  conclude  that  our  model  works 
for  samples  levels  of  roughness  for  1.5  nm  or  lower,  roughly  15%  of 
well-width.  It  is  our  theory  that  with  increasing  IFR  the  charge 
carrier  localization  builds  up  enough  so  that  a  breakdown  of  the 
quantum-well  effects  occurs,  which  results  in  an  exponential  in¬ 
crease  in  electrical  resistivity,  i.e.,  Sample  1019_c.  It  is  therefore  evi¬ 
dent  that  consistent  and  uniform  low  levels  of  roughness  within  the 
multilayer  QW  structure  are  of  high  concern  when  scaling-up  of 
this  technology.  The  roughness  values  shown  here  may  appear 
large  in  comparison  to  previously  reported  numbers  of  0.1  to 
0.4  nm  for  GaN/AlGaN  films  [14],  however,  this  is  attributed  to 
the  fact  that  these  films  were  grown  by  sputtering,  not  epitaxially. 
In  addition,  these  values  represent  the  average  z-roughness  over 
the  full  50  periods,  and  not  just  a  single  quantum  well. 

As  compared  to  commercially  available  n-type  PbTe-based  mate¬ 
rials,  which  have  a  thermoelectric  power  factor  between  0.0008  and 
0.0012  W/m2  K  over  the  temperature  range  of  500-700  K  [21  ],  our 
results  indicate  improved  properties  ranging  from  0.0018  to 
0.0023  over  the  same  temperature  range.  In  order  to  draw  further 
conclusions  on  material  performance  however  we  must  obtain  ther¬ 
mal  conductivity  measurements  at  high  operating  temperatures.  In 
a  separate  study  on  Si/SiC  multilayer  QW  structures  the  in-plane 
thermal  conductivity  showed  no  variation  from  that  of  the  bulk 
properties  for  temperatures  up  to  300  K  [22].  Based  on  this  assump¬ 
tion  that  our  thermal  conductivity  is  a  combination  of  the  parent 
bulk  material  properties  we  have  then  predicted  a  ZT  value  of 


Table  1 

An  XRR  analysis  summary  of  deposited  samples.  Measured  well-widths  are  within  2% 
of  the  target,  and  all  samples  have  low  IFR.  The  sample  with  10%  indium  doping 
represents  non-significant  data  due  to  the  lack  of  layer  periodicity  induced  from  a 
high  level  of  film  roughness. 


Sample 

Periods 

Width  (nm) 

IFR  (nm) 

AZO/InAZ02 

50 

9.8 

1.0 

AZO/InAZ05 

50 

10.2 

1.3 

AZO/InAZ08 

50 

10.1 

1.4 

AZO/InAZOlO 

50 

24 

23 
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Table  2 

Summary  of  IFR  values  for  four  50  period  AZO/Ino.osAZO  multilayer  structures  and  their  respective  measured  and  calculated  electrical  resistivity  at  325  K. 


Sample  ID 

Details 

W  (nm) 

A  (nm) 

A  (nm) 

A/w 

Pcaic.  (mQ-cm) 

Pmeas.  (mQ-Cm) 

%  Error 

1 009_a 

50  periods/AZO/In0.05AZO 

10. 

1.4 

7.8 

0.14 

10.4 

9.5 

8.6 

1 009_b 

50  periods/AZO/In0.05AZO 

9.9 

2.1 

19.5 

0.21 

44.4 

42.3 

4.7 

1019_c 

50  periods/AZO/In0.05AZO 

10.1 

4.1 

33.2 

0.41 

343 

413 

20.4 

1025_b 

50  periods/AZO/In0.05AZO 

10.2 

1.3 

8.9 

0.13 

10.2 

9.4 

7.8 

1 .3-1.5  for  the  AZO/Ino.osAZO  sample  at  600  K.  Other  studies  how¬ 
ever  have  shown  that  both  the  multilayer  interfaces  and  poly-crys¬ 
talline  grain  boundaries  will  act  as  scattering  mechanisms  for 
phonons  [23],  which  makes  it  possible  that  the  QW  structure  will 
reduce  the  thermal  conductivity,  which  would  ultimately  result  in 
even  a  larger  ZT. 

5.  Conclusions 

In  conclusion,  high  quality  multilayer  structures  of  Al-doped 
ZnO  and  (In,  Al)  co-doped  ZnO  layers  have  been  deposited  using 
RF  sputtering.  The  fabricated  multilayer  quantum  well  structures 
deposited  show  a  reduction  in  electrical  resistivity  by  an  order  of 
magnitude  at  2%  and  5%  indium  well  doping  as  compared  to  the 
bulk  films.  This  combined  with  material  dependent  Seebeck  coeffi¬ 
cient  enhancements  results  in  a  3x  increase  in  thermoelectric 
power  factor  and  an  estimated  ZT  of  1.3-1. 5  for  the  AZO/Ino.osAZO 
multilayer  quantum  well  sample  at  600  K.  In  addition,  the  local 
flatness  of  the  multilayer  structure  has  proven  to  be  critical  for 
achieving  optimum  thermoelectric  performance.  Multilayer  sam¬ 
ples  deposited  resulted  in  a  wide  range  of  IFR  ranging  from  10% 
to  40%  of  the  well-width.  From  experimental  and  theoretical  calcu¬ 
lations  it  is  desired  to  have  an  IFR  to  well-width  ratio  <0.1  in  order 
to  have  optimized  thermoelectric  transport  within  the  multilayer 
structure. 
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